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ABSTRACT 

The discovery of two accreting millisecond X-ray pulsars in binaries with « 43 minute orbital periods allows 
for a new probe of the donor's structure. For XTE J1751-305, only a hot white dwarf (WD) can fill the Roche 
Lobe. A cold He WD is a possible solution for XTE J0929-314, though I will show that evolutionary arguments 
make a hot WD more likely. In addition to being larger than the T = models, these finite entropy, low-mass 
(Mc < O.O3M0) WDs have a minimum mass for a fixed core temperature. If they remain hot as they lose mass 
and expand, they can "evaporate" to leave an isolated millisecond radio pulsar. They also adiabatically expand 
upon mass loss at a rate faster than the growth of the Roche radius if the angular momentum deposited in the disk 
is not returned to the donor. If the timescale of the resulting runaway mass transfer is shorter than the viscous 
timescale in the outer disk, then the mass transfer instability of Ruderman and Shaham for He WDs would be 
realized. However, my estimates of these timescales still makes the instability unlikely for adiabatic responses. I 
close by noting the possible impact of finite T WDs on our understanding of AM CVn binaries. 
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1. INTRODUCTION 

The recent discovery by the Rossi X-Ray Timing Explorer 
of two accreting millisecond pulsar transients XTE J 175 1-305 
{vs « 435 Hz; Markwardt et al. 2002) and XTE J0929-314 
{vs « 185 Hz; Galloway et al. 2002) has allowed us to learn 
about the donors in these ultracompact binaries. At forb ~ 43 
minute orbital periods, H-rich donors are ruled out (Nelson, 
Rappaport & Joss 1986), and He-rich stars (see Podsiadlowski, 
Rappaport & Pfahl 2002 for an updated discussion of the range 
of H/He ratios) or WDs are filling the Roche lobe (RL). If a de- 
generate WD, then the orbital period increases as mass is trans- 
ferred at the rate set by angular momentum losses from gravity 
waves, 7gr (see Verbunt 1993). 

The measured pulsar orbital parameters yield the RL filling 
companion's mass, Mc, and the minimum mass transfer rate, 
Mgr = 3McJcr/2J, which are shown in Figure 1 for neutron 
stars of = 1.4-2.OM0. The low M's are the likely cause 
for the transient behavior, as the steady-state outer disk temper- 
ature is below these element's ionization temperature (Tsug- 
awa & Osaki 1997; Menou et al. 2002), even if X-ray heating 
is included at the rate inferred in other X-ray binaries (Dubus, 
Hameury & Lasota 2001). If identical WDs are the donors in 
both of these binaries, then their identical orbital periods would 
require the same RL filling solutions. This would constrain 
the inclination for XTE J0929-314 to less than 37 degrees and 
M, > O.OBMq for both XTE J0929-314 and XTE J1751-305. 

The dotted lines in Figure 2 show the mass-radius rela- 
tion for RL filling donors of XTE J 175 1-305 (lower line) and 
XTE J0929-3 14 (upper line). The solid and dashed lines (which 
exhibit a maximum radius due to the onset of Coulomb physics) 
are the cold (r = 0) WDs of Zapolsky & Salpeter (1969, here- 
after ZS) for pure He and C, respectively. A cold He WD will 
fill the RL for XTE J0929-3 14 (Galloway et al. 2002), whereas 
there are no cold WD solutions that fill the RL for XTE J 1751- 
305 (Markwardt et al. 2002). 



This evidence for finite 7^. WDs motivated my calculations 
shown in Figure 2 by the solid and dashed lines that diverge 
at low Mc- These models (see §2) are for He at = 10^ K 
and 10^ K and C at 7; = lO*" K and T, = 3 x 10^ K. These 
provide RL filling solutions, but don't a priori differentiate be- 
tween He or C or any other element that might dominate the 
donor star. I consider C WDs throughout this paper, as Schulz 
et al. (2001) have measured a high Ne to O ratio in the matter 
transferred onto NSs in ultracompact X-ray binaries (confirmed 
by Juett, Psaltis & Chakrabaity 2002; Homer et al. 2002; Juett 
& Chakrabarty 2002). These measurements led Schulz et al. 
(2001) and Juett et al. (2001) to suggest that the donors in these 
binaries are the cores of previously crystallized C/O WDs. 

In §2 1 justify (on evolutionary grounds) that hot WDs are ex- 
pected, describe the construction of their RdMc) relations, and 
point out the existence of a minimum Mc solution for a given 
Tc- I put the models in the context of mass donors in §3 and 
calculate 

d\nRc 
dmMc 

the adiabatic exponent needed to evaluate the stability of mass 
transfer. Though uncertainties remain regarding the ability 
of the accretion disk to store angular momentum on long 
timescales (Verbunt & Rappaport 1988), my initial work finds 
that the expansion of a He WD donor under mass transfer can 
exceed that of the RL, possibly allowing for the Ruderman & 
Shaham (1983) instabihty. 

2. FINITE ENTROPY WHITE DWARFS AND EVAPORATION 

The material that fills the RL in these systems today was 
once the deep interior of an initially more massive WD, so I 
begin by discussing the extent to which prior evolution (before 
RL filling) can affect the temperature of the current low-mass 
WDs. There are a few evolutionary scenarios that place a He or 
C/O WD in a tight enough orbit (following a common-envelope 
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induced spiral-in) about a NS that gravity wave emission will 
place it in contact within 5-10 Gyr (Rasio, Pfahl & Rappaport 
2000; Yungelson, Nelemans & van den Heuvel 2002). In this 
case the WD has had some time to cool before initiating RL 
overflow. Nelemans et al. (2001) discuss a similar scenario for 
the origin of AM CVn binaries. 

The rate at which they cool differentiates He and C/O WDs. 
The larger specific heat of a He WD slows its cooling (Althaus 
& Benvenuto 1997; Hansen & Phinney 1998). For example, a 
O.2OM0 He WD would have 7; = 10(3.3) x lO*" K at « 1.0(4.0) 
Gyr (Althaus & Benvenuto 1997) and be liquid when RL filling 
occurs. The smaller specific heat of a Q.6Mq C/O WD allows it 
to cool to 2.5 X lO^K in 4 Gyr and begin crystallization (Salaris 
et al. 2000). Hence, most C/O WDs would be solid prior to RL 
filUng (this is needed for the enhancement of ^^Ne by fraction- 
ation; Schulz et al. 2001) and might only melt if heated during 
the GW inspiral (Iben, Tutukov & Fedorova 1998). The other 
differentiation between such WDs comes from their mass, as 
Yungelson et al. (2002) argue that the mass transfer at the on- 
set of RL fining will be unstable if Mc > O.44M0, excluding 
traditional C/O WDs from ever reaching short orbital periods. 
However, Yungelson et al. (2002) note that "hybrid" WDs with 
C/O cores can have Mc < 0.44Mq and thus evolve to lower 
mass systems filling the RL at 20-40 minutes. 

The inability for the donor WDs to cool on the ^ Gyr 
timescale of the mass transfer phase (Rappaport et al. 1987) 
makes it clear that the initial entropy is the minimum value at- 
tainable. All of the entropy is in the liquid ions, which have 
an adiabatic scaling, Tc oc (Hernanz et al. 1988), less 

steep than for an ideal gas (T oc p^^^). Since pc oc M^, a 
« 0.013Mq He WD made by adiabatically expanding the deep 
interior of a 0.20Mq WD would be a factor of « 15 cooler, 
or Tc ~ 6.6(2.2) x lO^K for the initial range of T's discussed 
earlier Of course, any tidal heating will increase Tc (see §3). 

So I have constructed Mc < 0.03Mo WDs of finite Tc that 
contains the Coulomb physics (from Chabrier & Potekhin 
1988) that causes the turnover to constant density "rocks" in 
the mass-radius relation at « 1Q~^ Mq seen in the ZS models in 
Figure 2. I have imposed an arbitrary (but convectively stable) 
temperature profile ofT = TdP/Pcr^ while integrating hydro- 
static balance and mass conservation. I used the electron equa- 
tion of state of Paczynski (1983) and halt integrations at the 
point where the pressure has fallen to 10~^ of the central value, 
which avoids the need for an envelope model. The results are 
shown in Figure 2 and make clear that the temperatures ex- 
pected from adiabatic evolution of a He WD are adequate to 
provide the slight radius expansion needed to fill the RL for 
these accreting millisecond pulsars. This is easier for He WDs 
because of the extra number of ions per electron compared to 
C/O WDs. White dwarfs made of C/O will need some tidal 
heating to fill these RLs. 

I also found that there is a minimum mass solution for a WD 
of fixed Tc- This is related to the electron Fermi energy of the 
nearly constant density solutions at low masses and is reflected 
in the radius divergence of the models in Figure 2 at low Mc- 
My semi-analytic modeling finds that these finite Tc solutions 
eventually "turnover" to follow a track of Rc oc Mc, as expected 
for an ideal gas polytrope (where Tc oc Mc/Rc)- Since there 
is no way to reach these solutions of higher Mc while mass is 
being lost, these "evaporative" endpoints where dlnRc/dlnMc 
diverges can only be reached if tidal heating is adequate to keep 
these models hot while the central density is dropping. More 



work remains to actually show that this will eliminate the com- 
panion and lead to an isolated millisecond radio pulsar. 

3. MASS TRANSFER INSTABILITIES AND QUIESCENT EMISSION 

Barring tidal heating adequate to reach the "evaporative" 
endpoint alluded to in §2, I will now discuss a possible mass 
transfer instability. Matter leaving the WD takes angular mo- 
mentum with it as it settles into the accretion disk. This an- 
gular momentum is returned to the WD via tidal torques once 
there has been enough time for viscosity to move the material 
outwards from the splash point. The timescale for this angu- 
lar momentum loop to be closed has been considered to be long 
enough (or uncertain enough) that Ruderman & Shaham (1983) 
raised the possibility that a mass transfer instability could oc- 
cur when the expansion of a low-mass donor (Mc ^ Mj) due to 
mass loss (measured by n^d of eq. |[lj) exceeds that of the RL, 
measured by 

dlnRu d\na 1 

dlnMc " dlnMc^3' *• "* 

where a is the orbital separation. The value of mr presuming 
the angular momentum "sink" of the disk (using the fitting for- 
mula of Verbunt & Rappaport 1988) is shown in Figure 3 by the 
dotted lines for (from left to right) Mj = 0.6, 1 .0, 1 .4, & 1 .8 Mq . 

Ruderman & Shaham (1983) presumed n^d = -1/3, which 
was shown to be an overestimate by Hut & Paczynski (1984) 
and Bonsema & van den Heuvel (1985), who found n from the 
ZS relation (shown in Figure 3 by the solid[dashed] lines for 
He[C]). These authors noted that a ZS model moved the cross- 
ing (wAd < wr) to such low Mc that the binary might not reach 
it in a Hubble time. However, for finite T He WDs, I show 
that this conclusion is altered. The solid circles are my evalua- 
tion of Had for the = 10^ K (set closest to the solid line) and 
Tc = 10^ K (bottom-most set) He WDs of Figure 2 and show 
that the finite entropy model sits between the T = Q value and 
the "naive" n = -1/3 guess for a perfect gas. For a M^- = 1 .4 Mq 
accretor, the instability question is first raised (since n^d < mr) 
when a 10^ K He WD has Mc « 0.01 Mq. This will occur 
about a Gyr after mass transfer has started presuming just grav- 
ity wave emission. My evaluations of /lAd for C WDs at the r^'s 
of Figure 2 never found such a crossing. 

The instability question is thus raised for He WDs and re- 
quires some extra entropy from tidal heating. The amount 
of heating is small, as for liquid He to reach Tc ^ lO*" K 
from about 1/2 that value (which it had initially) requires 
3kB(0.5Tc)/4mp « 3 x 10'^ erg g"' (the specific heat in the hq- 
uid state is « 3kBT per ion). This energy is comparable to the 
current rotational energy per gram for a tidally locked WD at 
the 40 minute orbital period. Hence, only a fraction of the 
higher rotational energy per gram from the tidally locked WD at 
a shorter orbital period needs to be deposited to yield « 10^ K 
today. 

The remaining criticism against the mass transfer instability 
is one of relative timescales (Verbunt & Rappaport 1988; Pried- 
horsky & Verbunt 1988). The mass-transfer instability most 
likely grows on a timescale Tg sa McH / McRc (Verbunt & Rappa- 
port 1988), where H is the scale height in the WD atmosphere, 
fixed by the X-ray heating in quiescence from the hot NS (Bild- 
sten & Chakrabarty 2001). The binary spends most of the time 
undergoing mass transfer to the outer disk, but little to no ac- 
cretion onto the NS. Bildsten & Chakrabarty (2001) showed 
that the X-ray emission always detected from NS binaries in 
quiescence (see Bildsten & Rutledge 2001 for an overview) 
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heat the companion on the side facing the NS. Presuming the 
NS thermal emission is at the level predicted by Brown, Bild- 
sten & Rutledge (1998) (giving Lx.q « 5 x lO"'^ erg s"') then 
the WDs in XTE J1751-305 and XTE J0929-314 will have 
Tm ~ 5000-6500 K.' The scale height is H/Rc « 3 x 
yielding Tj « 10* yr. This is still much longer than the time 
for viscosity to finally play a role in moving material outwards 
in the outer accretion disk, most likely allowing for the angular 
momentum to get back to the donor and move it out, so that 
hr = -5/3. The adiabatic index is never less than -5/3. How- 
ever, the isothermal WD response is < -5/3 and mass transfer 
would be unstable for Mc < 4(9.5) x \Q~^Mq in a He WD of 
Tc = 4(10) X 10^ K as long as internal heating can keep the WD 
isothermal as mass is lost. These masses are near the values 
I noted earlier for the "evaporation", namely where the mass- 
radius relation turns over 

4. CONCLUSIONS 

The timing of two accreting millisecond pulsars (Markwardt 
et al. 2002; Galloway et al. 2002) in ultracompact binaries has 
probed the WD donor properties to new levels and shown that 
they are of finite entropy. This motivated my calculations of 
low-mass WDs of finite Tc that allow for Tc to be constrained. 
For He WDs, the implied Tc are nearly that expected just from 
adiabatic expansion of the initially hot WD that filled the RL. 
Only a small amount of tidal heating is needed. More tidal 
heating is needed to make a C/O WD fill the RL. 

These finite T solutions allowed for a re-evaluation of Ruder- 
man & Shaham's (1983) scenario for making isolated millisec- 



ond radio pulsars via a mass transfer instability. I find that the 
adiabatic mass transfer instability can occur for a hot He WD 
as long as the angular momentum leaving the RL filling star is 
not returned to the orbit. I have thus eliminated one criticism of 
their model, though the question of angular momentum elim- 
ination remains a serious one. The physics of hot, low-mass 
WDs also yields a minimum mass WD solution for a fixed Tc 
so that a mass transfer instability can occur if the donor remains 
isothermal under mass loss. I thus speculate that an evaporative 
or mass transfer instability endpoint might occur as long as tidal 
(or other) heating persists at 40-80 minute orbital periods. 

This work also impacts AM CVn binaries, where a low mass 
He star donates material to a more massive WD (see Solheim 
1995 for a review). The larger WD radii lead to more GW emis- 
sion and a higher Mc than expected for a given Porb- Hence, 
models which track the WD entropy will fall between the de- 
generate and non-degenerate models in Nelemans et al. (2001) 
and depend on both the age of the system when RL filling oc- 
curs (as this fixes the initial WD entropy) and any tidal heating 
that occurs during the mass transfer If either of the instabilities 
discussed above occur, then the endpoint of AM CVn's could 
well be a DB WD (e.g. Tutukov & Yungelson 1996). 

I thank Deepto Chakrabarty for alerting me to the discovery 
of these transients and for many conversations. Ira Wasserman 
provided great physics insights during the progress of this work, 
which was supported by NASA via grant NAG 5-8658 and by 
the NSF under grants PHY99-07949 and ASTOl-9642. L. B. is 
a Cottrell Scholar of the Research Corporation. 
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'These high TEff's should aid detection of the quiescent counterpart, as the luminosity (presuming Rc = 0.045Rq) will be Lc ~ (0.5 — 1.0) X 10 ^Lq just from 
reprocessing the NS's thermal emission. Presuming rgff « 6000 K and using the WD radius implied by RL filling, I get Mv 12 for the heated WD face. For the 
high latitude source XTE J0929-314 at d Ri 5-7 kpc (Galloway et al. 2002), this gives V 26. 
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Fig. 1 . — The companion masses and gravity-wave driven mass transfer rates for XTE J1751-305 and XTE J0929-314 as a function of cosine of the inclination 
angle. The hatched regions (upper is XTE J1751-305 and lower is XTE J0929-314) are bounded below by = 1.4Mq and above hy = 2.0Mq. It was presumed 
that the donor responds to mass transfer with n = -1/3, only a small difference from riAi- 
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Fig. 2.— Roche lobe filling mass-radius relations for XTE J1751-305 and XTE J0929-314 compared to WDs. The dotted lines denote the RL filling donors 
of XTE J1751-305 (lower line) and XTE J0929-314 (upper line). The soUd (dashed) lines are T = WDs of ZS for pure He (C) using Rappaport & Joss's (1984) 
correction to ZS. The lines which diverge at small Mc are hot WDs. The He WDs have Tc = K and 10* K and the C WDs have 7^ = 10* K and 7^ = 3 x 10* K. 
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Fig. 3. — Roche radius and stellar radius changes with mass transfer. The dotted lines show nj? for (from left to right) Mx = 0.6, 1.0, 1.4 & I.SMq when the 
transferred angular momentum is stored in the disk. The solid (dashed) line are WAd = "zs for 7 = pure He (C) WDs. The filled circles are nAd for He WDs of 
Ti = 10^ K (set closest to the soUd line) and Tc = 10* K. 



